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A search for leptoquarks has been performed in 310 pb -1 of data from pp collisions at a center- 
of-mass energy of 1.96 TeV, collected by the DO detector at the Fermilab Tevatron Collider. The 



4 



topology analyzed consists of acoplanar jets with missing transverse energy. The data show good 
agreement with standard model expectations, and a lower mass limit of 136 GeV has been set at 
the 95% C.L. for a scalar leptoquark decaying exclusively into a quark and a neutrino. 

PACS numbers: 14.80.-j, 13.85.Rm 



Many extensions of the standard model (SM) that at- 
tempt to explain the apparent symmetry between quarks 
and leptons predict the existence of leptoquarks (LQ) . 
These new particles are scalar or vector bosons that carry 
the quantum numbers of a quark-lepton system. They 
are expected to decay into a quark and a charged lep- 
ton with a branching fraction (3, or into a quark and a 
neutrino with a branching fraction (1—0). At pp collid- 
ers, leptoquarks can be pair produced, if sufficiently light, 
primarily by qq annihilation and gluon-gluon fusion, with 
a production cross section independent of the unknown 
leptoquark-quark-lepton coupling. For (3 = 0, the result- 
ing final state consists of a pair of acoplanar quark jets 
with missing transverse energy, 1£t, carried away by the 
two neutrinos. 

In this Letter, a search for leptoquarks that decay into 
a quark and a neutrino, using data collected at a center- 
of-mass energy of 1.96 TeV with the DO detector during 
Run II of the Fermilab Tevatron Collider, is reported. 
The production cross section for vector leptoquark pairs 
is larger than that for scalar leptoquarks, but it is model 
dependent. The interpretation of the results is therefore 
presented in terms of scalar leptoquark masses. The most 
constraining 95% C.L. lower mass limit for (3 = 0, pre- 
vious to this search, was 117 GeV, obtained by the CDF 
Collaboration with 191pb _1 of Run II data|2|]. 

A detailed description of the DO detector can be found 
in Ref. Q. The central tracking system consists of a 
silicon microstrip tracker and a fiber tracker, both lo- 
cated within a 2 T superconducting solenoidal magnet. 
A liquid-argon and uranium calorimeter covers pseudora- 
pidities up to |^| w 4.2, where r\ = — In [tan (0/2)] and 9 
is the polar angle with respect to the proton beam direc- 
tion. The calorimeter consists of three sections housed 
in separate cryostats: the central one covers \r]\ < 1.1, 
and the two end sections extend the coverage to larger 
The calorimeter is segmented in depth, with four 
electromagnetic layers followed by up to five hadronic 
layers. It is also segmented in projective towers of size 
0.1 x 0.1 in r/-(j) space, where (f> is the azimuthal angle 
in radians. Calorimeter cells are formed by the intersec- 
tions of towers and layers. Additional energy sampling 
is provided by scintillating tiles between cryostats. An 
outer muon system, covering \r]\ < 2, consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar layers 
beyond the toroids. 

For this search, data collected with a jets + Jpr trigger 
have been analyzed. At the first level, this trigger selects 
events in which at least three calorimeter trigger towers 



of size A(f> x A?7 = 0.2 x 0.2 record a transverse energy 
in excess of 5 GeV. At the second and third trigger lev- 
els, requirements are placed oiiI/It, the vector sum of the 
jet transverse momenta (I/It — | Sjcts'PVI)- Coarse jets 
are reconstructed from trigger towers at the second level, 
while the full detector information is used at the third 
level. The I/It thresholds are 20 and 30 GeV at the sec- 
ond and third levels, respectively. The trigger efficiency 
is larger than 98% for events fulfilling the selection cri- 
teria of this analysis. Data quality requirements on the 
performance of each detector subsystem yielded an avail- 
able integrated luminosity of 310 pb -1 . 

The offline analysis utilized jets reconstructed with 
the iterative midpoint cone algorithm with a cone 
size of 0.5. The jet energy calibration was derived from 
the transverse momentum balance in photon+jet events. 
Only jets with pt > 15 GeV that passed general qual- 
ity criteria, based on the jet longitudinal profile in the 
calorimeter, were selected for this analysis. The mis- 
sing transverse energy was calculated from all calorimeter 
cells, corrected for the energy calibration of reconstructed 
jets and for the momentum of reconstructed muons. 

The sample of approximately 14 million events col- 
lected with the jets +$t trigger was reduced by requiring 
the following preselection criteria to be satisfied: at least 
two jets; $t > 40 GeV; If T > 40 GeV, where, in con- 
trast to I/Iti information from energy not belonging to 
reconstructed jets is taken into account; and A$ < 165°, 
where A$ is the acoplanarity of the two leading jets, i.e., 
the two jets with the largest transverse momenta, defined 
as the difference between their azimuthal angles. To en- 
sure that the selected events were well contained in the 
detector, the position of the interaction vertex along the 
beam direction was required to be within 60 cm of the 
detector center. 

Events in which the presence of obvious calorimeter 
noise could be detected were rejected. The inefficiency as- 
sociated with this procedure was measured using events 
collected at random beam crossings (zero- bias events), 
and events collected with an unbiased trigger and con- 
taining exactly two jets back-to-back in azimuth. At this 
stage, 306,937 events survived. 

Signal efficiencies and SM backgrounds have been eval- 
uated using a full GEANT || based simulation of events, 
with a Poisson average of 0.8 minimum-bias events su- 
perimposed, corresponding to the luminosity profile in 
the data sample analyzed. These simulated events were 
reconstructed in the same way as the data. The jet 
energies further received calibration corrections and an 
additional smearing to take into account residual differ- 
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ences between data and simulation, as determined with 
photon+jet events. The instrumental background due to 
jet energy mismeasurements in QCD multijet production 
was estimated directly from the data. 

The SM processes expected to yield the largest back- 
ground contributions are vector boson production in as- 
sociation with jets, among which Z — ► vv is irreducible. 
Vector boson pair production and top quark production 
have also been considered. All of these processes were 
generated with alpgen 1.30, interfaced with pythia 
6.202 for the simulation of initial and final state ra- 
diation and for jet hadronization. The parton distribu- 
tion functions (PDFs) used were CTEQ5L 8]. The next- 
to-leading order (NLO) cross sections for vector boson 
production in association with jets were calculated with 
mcfm 3.4.4 9] and the CTEQ5M PDFs. 

The production of scalar leptoquarks via the processes 
qq or gg — > LQLQ was simulated with PYTHIA and the 
CTEQ5L PDFs. The chosen leptoquark masses ranged 
from 80 to 140 GeV, in steps of 5 GeV. For each mass, 
10,000 events were generated. The NLO leptoquark 
pair production cross sections were calculated using a 
program based on Ref. [l(J, with CTEQ6.1M PDFs [H|. 
For the mass range considered, they vary from 52.4 to 
2.38 pb. These nominal values were obtained for a renor- 
malization and factorization scale equal to the leptoquark 
mass. 

The selection criteria for this analysis are listed in Ta- 
ble [I] together with the numbers of events surviving at 
each step and with the cumulative efficiency for a lepto- 
quark mass of 140 GeV. The jet kinematic cuts CI to C4 
reject a large fraction of the SM and instrumental back- 
grounds. They take advantage of the central signal pro- 
duction and decay by requiring that |77dct | be smaller than 
1.5 for the two leading jets, where r]d e t is the pseudora- 
pidity measured from the detector center. Cut C5, where 
EMF is the fraction of jet energy contained in the elec- 
tromagnetic section of the calorimeter, rejects jets likely 
due to photons or electrons. 

In cut C6, the total transverse energy of the charged 
particles emanating from the interaction vertex and as- 
sociated with a jet, as measured in the tracking system, 
is compared to the jet transverse energy recorded in the 
calorimeter. The charged particle fraction CPF, i.e. the 
ratio of these two quantities, is expected to be close to 
zero either if a wrong interaction vertex was selected, in 
which case it is unlikely that the charged tracks truly as- 
sociated with the jet will come from the selected vertex, 
or if the jet is a fake one, e.g. due to calorimeter noise, 
in which case there should be no real charged tracks as- 
sociated with it. The efficiency of this jet confirmation 
procedure was determined using events containing two 
jets back-to-back in azimuth. 

Cut C7 was applied to suppress further the instrumen- 
tal background, which is enriched in multijet events by 
the acoplanarity requirement. The efficiency of such a jet 



TABLE I: Numbers of data events selected and signal cumu- 
lative efficiencies for ttilq — 140 GeV at various stages of the 
analysis. The leading and subleading jets are denoted jet-1 
and jet- 2. 



Cut applied 


Events left 


Signal eff. (%) 


Initial cuts 


306,937 


58.8 


CI: jet-1 p T > 60 GeV 


206,116 


48.7 


C2: jet-1 [rjdetl < 1.5 


160,323 


46.8 


C3: jet-2 p T > 50 GeV 


48,979 


24.8 


C4: jet-2 \ Vdct \ < 1.5 


42,028 


22.7 


C5: jet-1 jet-2 EMF < 0.95 


40,821 


22.3 


C6: jet-1 jet-2 CPF > 0.05 


34,746 


22.2 


C7: exactly two jets 


5,213 


15.3 


C8: $ T > 70 GeV 


492 


11.8 


C9: isolated electron veto 


465 


11.7 


CIO: isolated muon veto 


399 


11.6 


Cll: isolated track veto 


287 


10.0 


C12: A$ max - A$ min < 120° 


180 


9.4 


C13: A$ max + A$ min < 280° 


124 


8.4 


C14: JpT > 80 GeV 


86 


7.0 



multiplicity cut is sensitive to the modeling of initial and 
final state radiation (ISR/FSR). To verify the simulation 
of these effects, (Z — ► ee)+ > 2-jet events were selected 
in the data, and compared to a simulation by alpgen for 
the production of (Z — > ee)+2-jcts, with ISR/FSR jets 
added by pythia. The two leading jets were required 
to fulfil criteria similar to those used in the analysis, and 
the numbers of events with additional jets were compared 
between data and simulation, as well as the pr spectra of 
those jets. The small deficit observed in the simulation, 
located mostly at pr < 20 GeV, was used to correct the 
signal and background simulations, and the statistical 
power of this test was taken as a systematic uncertainty. 

After cut C8, the level of the instrumental background 
is largely reduced and is similar to the level of the SM 
backgrounds. The finally cut value (cut C14) was op- 
timized as explained below. 

Cuts C9, CIO and Cll, reject a large fraction of the 
events originating from W/Z+jet processes. In cut C9, 
an electron with pt > 10 GeV is declared isolated if the 
calorimeter energy in a cone of radius 0.4 in ry-cj) around 
the electron direction does not exceed the energy con- 
tained in the electromagnetic layers inside a cone of ra- 
dius 0.2 by more than 15%. In cut CIO, a muon with 
Pt > 10 GeV is declared isolated if the calorimeter energy 
in a hollow cone with inner and outer radii 0.1 and 0.4 
around the muon direction is smaller than 2.5 GeV, and if 
the sum of the transverse energies of charged tracks, other 
than the muon, in a cone of radius 0.5 is smaller than 
2.5 GeV. In cut Cll, a charged track with px > 5 GeV is 
declared isolated if no charged track with pt > 0.5 GeV 
is found within a hollow cone of radii 0.1 and 0.4 around 
the track considered. This cut was specifically designed 
to reject hadronic decays of r-leptons; the use of a hol- 
low, rather than full cone renders it efficient also in case 
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FIG. 1: Distributions of A<E> max — A<E> min (a) and of A$ max + 
A$ m in (b) for data (points with error bars), for SM back- 
grounds (shaded histograms) , and for a 140 GeV LQ sig- 
nal (hatched histograms). In the A<I> max — A$ m i n distribu- 
tion, cuts CI to Cll are applied. The excess in data be- 
yond 120° is attributed to the non-simulated instrumental 
background. In the A$ max + A$ m i n distribution, the cut 
A$ max — A<J? m i n < 120° (C12) has been applied in addition. 
The locations of cuts C12 and C13 are indicated by arrows 
in (a) and (b), respectively. 



FIG. 2: Distributions of IpT for data (points with error 
bars), for SM backgrounds (heavy-shaded histograms), for 
the instrumental background (labeled QCD, light-shaded his- 
tograms), and for a 140 GeV LQ signal. In (a), all cuts except 
C8 and C14 are applied, the LQ signal is shown as a hatched 
histogram, and the insert shows how the instrumental back- 
ground contribution is estimated from power law (solid curve) 
and exponential (dashed curve) fits. The Ifir distribution in 
(b) is after all cuts, with the same shading code but with the 
signal contribution now displayed on top of all backgrounds. 



of decays into three charged particles. 

The angular correlations between the jet and IpT 
directions are used to suppress both the instrumen- 
tal and SM backgrounds. To this end, the minimum 
A$ m in($T, any jet) and maximum A$ max ($r, any jet) of 
the azimuthal angle differences between the direction 
and the direction of any of the two jets are combined as 
shown in Fig.^ It can be seen in Fig.QJi that cut C12 
rejects most of the remaining instrumental background, 
which is responsible for the excess beyond 120°. Cut 
C13, which suppresses SM backgrounds at the expense 
of a moderate reduction of the signal efficiency, was opti- 
mized as explained below. The variable A<i> max + A$ m ; n 
is the one which discriminates best the signal and the 
irreducible background from (Z — > w)+2-jets. Its effect 
is demonstrated in Fig.^i. 

Finally, the IpT and A<i> max + A$ min cuts were opti- 
mized for a 140 GeV LQ mass so as to minimize the cross 
section expected to be excluded in the absence of signal. 
Cut C8 was removed, and J^t cut values ranging from 



60 to 90 GeV were probed in 10 GeV steps. The cut on 
A$ max + A$ m j n was varied between 260° and 300° in 
steps of 10°. For each set of cuts, the instrumental back- 
ground was estimated as explained below. The system- 
atic uncertainties discussed further down were taken into 
account in the calculation of the expected limits. The 
optimal set of cuts reported as C13 and C14 in Table |U 
selects 86 data events. 

The instrumental background was estimated from ex- 
ponential and power law fits to the $t distribution (in- 
sert of Fig.EK) in the range [40, 60] GeV, where the sig- 
nal contribution is negligible, after subtraction of the 
SM expectation. Both fits were extrapolated beyond 
the JfiT cut value, and the average of the two results was 
taken as the instrumental background estimate, with a 
systematic uncertainty accounting for the difference be- 
tween the two fit results. The final Ipx distribution is 
shown in Fig.|2t>. The values of the SM and instrumental 
backgrounds are given in Table |nj The largest back- 
ground sources are, as expected, (Z — > vv)+2-]ets and 
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TABLE II: Numbers of events expected from standard model, 
instrumental and total backgrounds; number of data events 
selected; and number of signal events expected for tulq ~ 
140 GeV, assuming the nominal production cross section. For 
the total SM and total backgrounds, as well as for the signal, 
the first uncertainties are statistical and the second system- 
atic. The uncertainties on the individual SM backgrounds are 
statistical. The uncertainty on the instrumental background 
is mostly systematic from the difference between the power 
law and exponential fits. 



(Z -> vv)+2-jets 
(W -> fe)±jets 
(Z -> «)+jcts 
ti 

WW, WZ, ZZ 



34.6 ±4.3 
35.01^ 
0.3j" 
1.9 ±0.1 
1.2 ±0.2 



3+0.4 
-0.2 



Total SM background 


79 Q +10.1 +10.6 
' z - 9 -9.7 -12.1 


Instrumental background 


2.3 ± 1.2 


Total background 


7 r 9 +10.1 +10.7 
-9.7 -12.2 


Data events selected 


86 


Signal (m LQ = 140 GeV) 


51.8 ± 1.8±|;| 



(W -> &/)+jets {£ = e, (j,, r). 

The signal efficiencies at various stages of the analysis 
are given in Table QJfor tulq — 140 GeV. The efficiency 
decreases together with the leptoquark mass, reaching 
1.6% at 100 GeV. The number of signal events expected 
for a leptoquark mass of 140 GeV is indicated in Table ITTI 

The following sources of systematic uncertainty are 
fully correlated between SM background and signal ex- 
pectations: the relative jet energy calibration between 
data and simulation: +g% for the SM background and 
_4% for the signal; the relative jet energy resolution be- 
tween data and simulation: ~^_\% for the SM background 
and negligible for the signal; the efficiency of the jet mul- 
tiplicity cut: ±3%, after corrections of —3% for the SM 
background and —2% for the signal; the trigger efficiency: 
±2% after all selection cuts; and the integrated luminos- 
ity of the analysis sample: ±6.5%. 

In addition to the 1]^% statistical uncertainty of the 
simulation, the normalization of the SM background ex- 
pectation is affected by a ±12% uncertainty, as inferred 
from a comparison of data and simulated (Z — > ee)±2-jet 
events selected with the same criteria for the jets as in the 
analysis sample. The uncertainty of ±1.2 events on the 
instrumental background was estimated from power law 
and exponential fits to the^r distribution, as explained 
previously. As a check, the same procedure was applied 
to the events with A$ max — A$ m j n > 120°, which are 
dominated by the instrumental background contribution. 
This showed that the high ]£t tail is somewhat underesti- 
mated, possibly by as much as nine events, which leads to 
conservative results in terms of limit setting. Finally, the 
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FIG. 3: Observed (solid curve) and expected (dashed curve) 
95% C.L. upper limits on the cross section times (1 — (3) 2 as 
functions of the leptoquark mass. The nominal cross section 
for scalar-leptoquark pair production is also shown for f3 = 
(dash-dotted curve), with the shaded band indicating the 
uncertainty due to the choices of PDFs and of renormalization 
and factorization scale. 



uncertainty on the signal efficiency due to the PDF choice 
was determined to be t^%, using the twenty-eigenvector 
basis of the CTEQ6.1M PDF set [n|. 

As can be seen in Table [H] and Fig.|2b, no signifi- 
cant excess of events is observed in the data above the 
background expectation. Therefore, given the number 
of selected events, the SM and instrumental background 
expectations, the integrated luminosity of 310pb _1 , the 
signal selection efficiency as a function of the leptoquark 
mass, and the statistical and systematic uncertainties dis- 
cussed above, a 95% C.L. upper limit on the cross section 
times (1 — (3) 2 has been determined as shown in Fig. [31 
using the modified frequentist C'L S approach 12]. The 
expected limit in the absence of signal is also indicated. 

The nominal theoretical cross section for the pair pro- 
duction of scalar leptoquarks is also shown in Fig.|3| It 
was obtained based on Ref. with CTEQ6.1M PDFs 
and for a renormalization and factorization scale /j, r f 
equal to the leptoquark mass. The uncertainty associ- 
ated with the PDF choice was estimated using the full 
set of CTEQ6.1M eigenvectors and combined quadrati- 
cally with the variations obtained when was modified 
by a factor of two up or down. For a leptoquark mass 
of 140 GeV, the PDF uncertainty on the theoretical cross 
section amounts to and the scale variation results 

in a change of the quadratic sum being ±%%. Re- 

ducing the nominal cross section by this theoretical un- 
certainty, shown as the shaded band in Fig-El a lower 
mass limit of 136 GeV is derived at the 95% C.L. Masses 
smaller than 85 GeV, to which this analysis is not sen- 
sitive, have been excluded previously 0, El • The cross 
section limit obtained here was combined with the results 
of the published DO search for first-generation scalar lep- 
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toquarks in the eeqq and evqq final states |14j , and the 
lower mass limit of 136 GeV was seen to hold independent 
of p. 

In summary, a search for acoplanar jet final states in 
pp collisions at 1.96 TeV, performed using a data sam- 
ple of 310 pb -1 collected by the DO detector, revealed no 
deviation from the standard model expectation. For a 
single-generation scalar leptoquark, a lower mass limit of 
136 GeV has been obtained for (3 = 0. While a tighter 
limit is available for third-generation leptoquarks |l5l | , 
due to the increased signal purity achieved with heavy 
flavor tagging, this is the most stringent limit to date for 
first- and second-generation scalar leptoquarks decaying 
exclusively into a quark and a neutrino. 
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